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Abstract 

We analyze azimuthal anisotropy in heavy ion collisions related to the reaction 
plane in terms of standard reggeon approach and find that it is nonzero even when 
the final state interaction is switched off. This effect can be interpreted in terms of 
Q ■ partonic structure of colliding nuclei. We use Feynman diagram analysis to describe 

. details of this mechanism. Main qualitative features of the appropriate azimuthal 

correlations are discussed. 

> ' 

in . 

^ Introduction 
p 

G^ I Investigations of azimuthal anisotropy of secondary particle distributions in heavy ion col- 
^ ■ lisions is used as an effective tool for analysis of dynamical mechanisms of these processes. 
O ■ Anisotropic flows and, in particular, the elliptic flow t>2 = (cos 20) are considered usually 
as an important source of information about properties of a dense state of hadronic matter 
\ formed in the anisotropic overlap region of colliding nuclei ^. Crucial point is the relation 
I of the anisotropy to the reaction plane orientation (azimuthal correlations not related to 
this orientation are called non-flow correlations). The dependencies of flow coefficients on 
c.o.m. energy, centrality, rapidity and transverse momentum of secondaries are ascribed 
usually to anisotropic properties of strongly interacting matter at early times after collision. 
It is assumed that anisotropic flows are due to final state interactions and in hydrodynamic 
models originate from pressure gradients connected to original asymmetry in the configura- 
tion space for noncentral collisions. Practically all existing models for anisotropic flows use 
classical description of processes. 

In this paper we use quantum relativistic approach, based on analysis of Feynman di- 
agrams and reggeon diagrams for a study of anisotropic flows in collisions of hadrons and 
nuclei at high energies. We show that a part of azimuthal anisotropy can be related to prop- 
erties of initial state of a fast nucleus. Namely, we analyze azimuthal correlations present in 
partonic configurations of fast colliding nuclei. Thus a part of the observed azimuthal asym- 
metry of final particles can come from parton correlations in initial state. We demonstrate 
how anisotropic flows in this case are related to the anisotropy of the overlap region. 

iSee e.g. [1] - [8]. 
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It is worth to point out that, though it is not possible to measure in general both 
transverse momenta and transverse coordinates of partons, but for fast moving heavy nucleus 
impact parameters of nucleons and fast partons can be fixed with a good accuracy. Therefore 
for colliding nuclei the radial directions in every point of transverse plane with respect to 
c.o.m. are well defined. Distributions of partons at each point are anisotropic in general. The 
collision process (if it is not absolutely central) selects partonic configurations asymmetric 
in impact parameter plane. This transforms the initial 'radial' asymmetry in colliding nuclei 
to the azimuthal asymmetry of produced particles momenta with respect to the collision 
plane. 

The simplest way to study such phenomena theoretically is to use the reggeon diagrams 
for inclusive cross sections, where all details of parton interactions are encoded in vertices 
and Regge trajectories - and we follow this way in the paper. 

The paper ^ is organized as follows. In Sec. 1 we calculate a contribution to V2 due to 
standard multiperipheral production described by reggeon diagram for inclusive production. 
The calculation is carried out in the momentum representation and gives f 2 as a function 
of impact parameter b and particle transverse momentum p^. A simple phenomenological 
parameterization is used for the vertex of a detected particle emission. Then in Sec. 2 we give 
a qualitative interpretation of the results in terms of the impact plane analysis. It is shown 
that the obtained azimuthal asymmetry can be ascribed to specific correlations presented 
at the level of partonic configuration of wave function of the fast nucleus. It exists without 
any final state interactions, which are usually considered as an obligatory condition for 
appearance of correlations related to the reaction plane orientation. We explain how these 
primary correlations manifest themselves as the elliptic flow due to anisotropy of the overlap 
region. In Sec. 3 we make an estimate of the anisotropy structure of the inclusive vertex in 
simple multiperipheral models and show that they naturally result in azimuthal dependence. 
Possible mechanisms enhancing the effect are discussed. In Sect. 4 we discuss generation of 
the directed flow effect in the same model. This effect is produced by nondiagonal in 
reggeons RP vertex and therefore it is concentrated at edges of rapidity space. In Sec. 5 we 
discuss contribution to 2-particle azimuthal correlations in reggeon approach and suggest 
possible tests of this mechanism. Sec. 6 contains discussion of multipomeron interactions 
which can have crucial influence on magnitude of the effect and on its atomic number and 
impact parameter dependencies. In Sec. 7 numerical estimates are presented to illustrate 
main qualitative features of the model. In Sec. 8 we briefly summarize our main results. 



1 Momentum representation 

Let us explain first the mechanism under consideration using reggeon diagrams. A simplest 
way to produce secondary particles at high energy is given by multiperipheral diagram 
of Fig. la. Inclusive cross section is described by the reggeon diagram of Fig. lb which 
corresponds to squared multiparticle amplitude with fixed momentum p of a single detected 
particle. 

Single ]:)articlc inchisivc cross section is related to the amplitude (j(q, p) of the diagram 

^ This paper is an extended version of a talk given at Session of Nuclear Physics Division of Russian 
Academy of Sciences in November 2007. 
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(a) (b) 
Figure 1: Reggeon diagram for one-particle inclusive production 



of Fig. lb at zero transfer momentum q = 0: 

Finci{pt;y) = ^^^^G{q^O;Pt,y), (1) 

where pt and y are transverse momentum and rapidity of the registered particle. We, 
however, need to calculate a contribution of this diagram at various values of q in order to 
find a correlation between transverse momentum pt and impact parameter value for colliding 
nuclei, b. This relation is given by Fourier transform of G(q; pt, y) 

FinciiPt, y; b) ~ G(b; p„ y) = j e^bq ^(q. ^) _ (2) 

Here 

^(q; Pi, y) = F^(g)L'(q; ?/)7(q, Pt)^(q; Y - y)FB{q) , (3) 

where Fa,b are 2-dimensional form factors of colliding nuclei which are Fourier transforms 
of their 2-dimensional profiles, TA^sip), 

FaMq) = / d^'b e-''^''TA,B{b) , ^^^(0) = A, Fb{0) = B, (4) 

and functions D {q; y) , D {q; Y — y) are connected to s-channel discontinuities of reggeon 
amplitudes ("cut reggeons") ^ 

L'(q; y) = 2 Im/(q; y) = ^^e-(^o+«'^)«^[«(o)-% . (5) 

An inclusive vertex of the reggeon diagram 7(q, p*) depends on 2-dimensional momenta pt 
and q. For estimates we will use the phenomenological parametrization 

7(q, Pt) = 7o exp{-rjq2 - r^p^ - er^iptq)^} . (6) 

Here dimensional coefficients r^, rp are expected to be of typical hadronic sizes about ro = 
1 fm, and e of order of 1. 

The dependence on variable (ptq) is crucial for appearance of azimuthal correlations 
tY]:)ical for elli])tic flow. Note that the terms proportional to the odd powers of p^q should 

■^We use an exponential parameterization of reggeon vertex g'jv exp(— iigq^) and linear Regge trajectory 
a(q2) = a(0) - a'q^. 



be absent if both reggeons have positive C-parity (hke pomerons) because the vertex should 
be even under the change pt ~Pt- More detailed analysis of structure and value of the 
vertex 7(q, p^) will be performed in sections 3,4. 

The elliptic flow f 2 is defined as a second Fourier coefficient in an azimuthal dependence 
of inclusive spectrum 

d 

ocl + ^ 2vn{b;pt, y) cos(n0) , (7) 

where is an azimuthal angle between impact parameter b which determines the reaction 
plane and transverse momentum p^. Thus, 

C2 ^ d(t) cos{2(j)) G{h]pt,y) 

V2{b;pt,y) = ^ r ^ ■ ^ 

Co J d(p G(b; pt,y) 

For the vertex function parameterized by eq.(6) one gets 

C^{b;pt,y) = j qdqUhq)FA{q^)FB{q^)D{q-y)D{q,Y -y)c^{q-h,p^^ n = 0,2, (9) 
with 

Co(g; h^pt) = H<pWm e-^^'^^'l-" , (10) 
C2(g;6,p,) = h{erlpiei2) e-o^'?'''/^ . 

where /q, i\ are the modified Bessel functions. 

To estimate an order of the effect let us first make the calculation with the approximate 
Gaussian parameterization of nuclear form factors FA^B(g^): 



Then 



FaM)=^M-^\b^'I^) ■ (11) 



-2(6;p*,z/) = ^i^ (12) 



with 



4 2/2 

tr^pib 
R^{R^ + 4e4pj) 



« = 02.02° 4 2, , R'^Rl + Rl + Rl,,m ■ (13) 



For a <^ 1 we have 



-i-(-)^(^)^(4)^- 

Thus, V2 is proportional to the "ellipticity" e of the inclusive vertex 7. The sign of V2 is 
the same as the one of e which can be both positive and negative depending on dynamical 
structure of the inclusive vertex (see discussion in Sect. 3). The formula (14) reveals an 
increase of V2 with transverse momentum p^ and impact parameter h. Results of numerical 
calculations are presented in sec. 7. 
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2 Impact parameter representation 



Presentation of formulas in terms of impact parameters allows us to give more instructive 
interpretation of the results. Eqs.(2),(3) in these variables have the form: 

^(b; pt, y)^j d'(3,j j d^h2 j d'^, Ta{(3,) h{h^ - f3,; y) 

X r(b2 - bi, Pi) - b2; r - y) TB{h - (3^) . (15) 

Here (see Fig.2) /3]^ and h — are the positions of participant nucleons of the nuclei A 
and B with nuclear profiles Ta and Tg, t(/3; y) is a Fourier transform of reggeon propagator 
D(q; y) and r(b; p*) is a Fourier transform of the inclusive vertex 7(q, pt): 



r(b;p,) 



|f,e-^(q,p.) 



(16) 



For the exponential model approximation (6), one has 
r(b;pt) = 



7re '^i'P* 



: exp 



62 



y 



(17) 



where is a component of b along p^ and by in the transverse direction. In order to reveal 





Figure 2: Interaction of two nucleons in impact plane 



a structure of Eq.(15) it can be written symbolically as a number of convolutions with 
notation f ® g ior convolution of functions / and g: 



FinciiPt, y, b) ~ Ta t{y) ®t{Y - y) ^Tb . 
Convolutions in eq.(15) can be done in arbitrary order, so it is equivalent to 

F..(P.y;b)^T..0r = /.2aT..(b-a)r(a;p.). 



(18) 



(19) 
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Here the function Tint which characterize a region of interaction of nuclei in b space is a 
convolution of the nuclear overlap * Tab = Ta®Tb and combined reggeon amplitude t regge- 

Tint ~ Tab ® tregge i 



where 



treggeih; Y) = t{y) ® t{Y - y) = ^^^f ^ ^""P ( 



regge\ / \ regge 



^Rle,,e{y) 



Rleggeiy) = '^Rl + C^'y , 



describes the regge-interaction amplitude of nucleons from colliding nuclei. 

For Gaussian parameterization of nuclear profiles Ta, Tb all convolutions have the Gaus- 
sian form again, resulting in the same final answer (12). 

In the impact parameter representation it is easy to see an origin of the correlation 
between directions of b and transverse momentum p^. Registration of the inclusive particle 
is realized by nonlocal probe which structure is determined by the inclusive vertex r(b; pt). 
This vertex has an eUiptic anisotropy in 6-plane along the direction of pt (see examples in 
Sec. 3). Due to this anisotropy the convolution of the overlap function T^nt and F is sensitive 
to gradients of nuclear densities. 

The size of the probe range is small in comparison to nuclear size Ra^Rb^ A^/^tq, so 
it is possible to calculate the integral (19) expanding the smooth distribution Ti„t(b — a) over 
small corrections a^/ RA^dy/ Ra- After integration over a^^ay we get explicitly anisotropic 
expression for the inclusive cross section (recall that the axis x is chosen along the momentum 
Pt). Omitting higher order terms in ai/RA we have in the approximation (6), (17) 

F^nei - TUh') + 2{2rl + ery,)^^^ + 4^^|^ [{rl + er^pX + r,6j] , (20) 

where derivatives are taken over the variable 6^. Thus, for the elliptic flow coefficient we 
have 

The second derivative of Tj„t(b^) contains a small parameter of order of (roZ-R^)^ which 
appears explicitly in Eq.(14). 

This formula can be interpreted in terms of parton probe. The local probe feels only 
parton density. An information about density gradients and its higher derivatives can be 
obtained by means of nonlocal probes with typically hadronic range. Being restricted only 
with lower multipoles a structure of the probe can be written as 

F(a; p) ~ (1 + hidi + + mpj)didj + . . . ) S^^^ (a) , (22) 

where di = d/dtti and 5 = Sr^, e = er^ are dimensional coefficients of typical hadronic range. 
The first term measures parton density as well as the isotropic one with 5ij, the second term 



* To avoid confusion let us stress that Tab describes how the number of nucleons from colliding nuclei 
which overlap in the impact plane depends on collision impact parameter b. It is an isotropic function not 
depending on the vector b direction. 
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feels gradients, and the e-term gives an information about ellipticity of partonic distribution 
at the point of production of registered particle. So the distribution over transverse momenta 
turns out to be correlated with impact parameter b: 

Finciiw, b) ~ Tint ® r = J d'^aTintih - a)r(a; 

^ TUb') + 25{p,h)TUb') + ^KVthfTiim , (23) 

where T' denotes the derivative over b^. This reproduces the elliptic flow coefficient V2 — 
e(ptb)^Tj'^^(6^) (cf. with (21)). The S-teim corresponds to the first flow coefficient Vi — 
5(ptb) T^^j(6^) (directed flow), which wc discuss in Sec. 4. 

Let us stress that partonic interpretation of reggeon diagrams depends on chosen Lorentz 
frame. Reggeon exchanges are highly non-local in space-time due to their complicated mul- 
tiparticle internal structure. From this viewpoint the fast hadron or nucleus is essentially 
a multiparton state of multiperipheral configuration (Fig. la represents the simplest exam- 
ple). The particles of such configuration are ordered in their rapidities and only most slow 
particles interact with a target. This means that in the lab frame related to the nucleus B 
the nucleus A is a multiparton state containing the detected particle with rapidity y and 
transverse momentum pt among others. This state has already an anisotropy in partonic 
distributions due to correlation between parton transverse momentum pt and its position 
bi in impact parameter plane. If nuclear collision is central then the overlap region is 
isotropic and after integration over bi anisotropy disappears. At non-central collision, how- 
ever, the overlap region has anisotropic almond shape and integration over this region keeps 
anisotropic distribution in p^. 



3 Inclusive vertex 

Let us estimate in simple models an anisotropy of inclusive pt probe related to the vertex 
7(q, Pt). If multiperipheral ladder of Fig.l consists of scalar particles with triple coupling 
then a structure of inclusive vertex is determined in lowest order by Fig. 3 where momentum 
Pt is fixed. Then the vertex can be written as 2-dimensional integral ^ 

7(q, Pt) ^g'j d'k G(k)G(k - pt)G(q - k)G(q - k + pt), (24) 

where particle propagators G depend on 2-dimensional momenta only. 

To estimate an anisotropy effect let us consider the vertex 7 at small q pt <^ jJ, (here 
~ Tq ^ is a typical mass parameter): 

7(pt,g,^)^7(p?)(l-e^^' + ---) , (25) 

A* 

^ There is an integration over 3-diniensional momentum k in the loop but for estimate we shall assume 

that longitudinal integration is already carried out. It introduces tmin terms into propagators which can be 
replaced by change of effective masses. So, in what follows we use propagators as functions of 2-dimensional 
momenta. 
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(a) (b) 

Figure 3: Simplest inclusive reggeon vertex 

where z — cos is a cosine of the angle between pt and q. The coefficient which determines 
a degree of ellipticity is denoted as e as we did before. For standard perturbative behavior 

of propagators: G'(k) = (k^ + numerical estimate gives a value e —0.2 for scalar 

couplings g in the diagram 3. Note that such a sign corresponds to negative (in contrast to 
experimental observation). This fact can be understood in a simple way by writing Eq.(24) 
in coordinate representation (see Fig. 3b): 

7(bi2, b34) I n ^^'b, G(6i3)G(632)G(624)G(64i) (26) 

i 

X 5'(bi2 - (b2 - bi))<52(b34 - (b4 - ba)) , 

where hij = bj — hj, and the vector bi2 is conjugated to q, while b34 - to the momentum 
Pf. The 2-dimensional Green functions G{h) ~ exp(— yu|b|) and therefore the integrand in 
(26) depends only on sum of all lengths |bi3| + |b32| + |b24| + |b4i|. It can be seen purely 
geometrically that maximum of ]^ C at fixed |b34|, |bi2| is reached in the configuration with 
parallel bi2 and b34. The maximal anisotropy (~ 0.3) corresponds to |bi2|/|b34| ~ 1 and 
decrease strongly with change of this ratio. Since 7(bi2, b34) is maximal for parallel bi2 and 
b34, then 7(q, p^) is maximal at parallel q and p^, i.e. e is negative. Similar arguments can 
be generalized to more complicated diagrams with spinless particles - configurations with 
parallel bi2 and b34 are more probable. 

The value and the sign of the effect depend on dynamical structure of the inclusive 
vertex 7(q, pt) - on form of propagators ^ and, more important, on structure of vertex of 
particle emission. Vertex properties are very sensitive to particle spin because it gives an 
extra azimuthal dependence. Calculations show that if the inclusive particle has spin 1 
(see, Fig. 4), then the sign of e is different (in accordance with data) and its value increases. 
Note that such situation is not far from reality because p-mesons give a considerable part 
of multipcriphcral particles and pions due to their decays have mainly close transverse 
momenta. Numerical estimate gives e ~ 0.5. 

Properties of the inclusive vertex 7(q, p^) can be studied experimentally in NN colli- 
sions, for instance by analysis of azimuthal two-particle correlations for hadrons in different 
windows of rapidity values. 

Let us discuss behaviour of the inchisivc vertex 7(q, p^) at large pt- In region p^ 3> /x^ 
it sliould decrease according to perturbation theory in power-like way, ~ (/^/pt)*", where 

^ Note that in case of Gaussian form of propagators, G{k) = exp(— fc^//i^), the product of propagators 
in coordinate space does not depend on the angle resulting in zero anisotropy, e = 0. 
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Figure 4: PP-diagram for inclusive p-meson production 



a degree m depends on dynamics. In simple scalar theory (24) 7(q, pt) ~ {f^/PtY, ^-nd 
2;-depending coefficients decrease even faster, leading to a small ellipticity: 

~ il^/Ptf at > 

It means that a size of the pj-probe becomes much smaller than /j,"^ and its form becomes 
more isotropic ^. 

The situation is different if the vertex includes particles with spin because in this case 
spin terms give an additional z dependence. As an example, in case of inclusive production 
of vector particles explicit calculation gives non- vanishing vertex ellipticity, e{pt) — > const 
when jj?. 



i{q ,pi,z = ^) -iir,pi,z = o) 



4 Directed flow 

The mechanism applied to generation of the elliptic flow V2 contributes to all i;„-harmonics. 
As it was mentioned in Sec.l there is a pecularity related to odd harmonics - the PP reggeon 
diagram do not contribute there because of symmetry of the PP inclusive vertex. The odd 
harmonics appear for vertices which couple different types of reggeons {RP vertices where 
P is a secondary reggeon, e.g. p) . 

Consider as an example contribution of diagrams of Fig. 5 to the first harmonic corre- 
sponding to the directed flow Vi. These diagrams contain the rapidity dependent factor 
which comes from the reggeon propagator Dji(q,y) ~ exp(Ap|/) where Ap ^ 0.5. The dia- 
grams (a) and (b) have different sign for the term qp^ so for the rapidity changing in the 
interval [— y, y] we get a factor fp{y) — 2exp(— A^y) sinh(Apy). 

We have to compare the term in proportional to qpt with the angle independent 

term 7(pt, q = 0) due to main diagram of Fig. lb (see Eq.(24)). For small pt 

l?''\ci, Pt)/7(Pt, q = 0) :^ AqP* ■ (27) 

It is clearly seen in the coordinate representation (26), since ~ ^34 ~ Pt"^ M^^j and the 
integrand depends on relative orientation of bi2 and b34 only in high orders in {n/pt) 
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Figure 5: Inclusive reggeon-pomeron diagrams contributing to vi. 



Calculation of the direct flow vi due to the RP diagrams are similar to ones made in Sect. 3 
and give (compare to (14)) 

v,c^S{roPt) g) (28) 

The numerical estimate of diagram of Fig. 6 gives 6 ~ 0.1 ^ 0.5. Qualitative behaviour of 
Vi as a function of 6 and pt is in agreement with experimental data (see Sect. 7). 




Figure 6: i?P-diagram for inclusive p-meson production 



Note also that in this model the expressions for high-order flows contain at small pt the 
same parameters as Vi and V2 but in higher powers. 

..-^»W(^)"(^)". (29) 

5 Two-particle azimuthal correlations 

We discussed hitherto correlations of transverse momentum of the detected particle with 
respect to the reaction plane. The reaction plane orientation is difficult to reconstruct from 
experimental data, and often another method for studying anisotropic flow is used. It is 
related with analysis of two-particle azimuthal distributions, which contain an information 
on anisotropic flow. As it is seen from analysis of experimental data the main part of 
these correlations is related to the Pt — b correlations for each particle. However, there 
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are possible sources of two-particle correlations unrelated to the reaction plane, so called 
'non-flow' correlations. 

Corresponding two-particle inclusive reggeon diagrams are shown in fig. 7. To calculate 
anisotropic flows these diagrams should be estimated not at zero momentum transfer q = 
as for standard bi-inclusive cross section but at fixed value of impact parameter b. 




(a) (b) (c) 



Figure 7: Reggeon diagrams for two particle inclusive production, (a) - particles p and k 
are emitted from different regions in transverse plane; (b) - p and k are emitted from a 
single reggeon chain; (c) - p and k are emitted from two reggeon chains attached to a single 
nucleon pair in colliding nuclei. 



The contribution of the diagram 7a is factorized in the impact parameter representation 
and contains factors corresponding to one-particle production of Fig. lb (see (23)). 



FlliPt, k,; b) ^ [TUb') + 4e(p, ■ hfTUb')] [T^nt{b') + • hrrUb') 
This results in factorized formula: 

(e2H0.,-0.,)) ^ J * - ' ^ V2{pt)v2{kt) . 

The contribution of diagram 7b is given by a convolution of Tint with two probes: 
F^^li{pt, kt; b) = y d^aid^a2 Tj„t(b - ai - a2)r(ai; Pt)r(a2; kj) . 



(30) 



(31) 



(32) 



This distribution (in the dipole approximation (22) with a — 0) will contain higher deriva- 
tives of Tint over 6^: 



F^^iiPt, kt; b) ~ T,nt + 86^(k, . pt)'T^, + 326^(k, ■ p,)(p, ■ b)(k, ■ b)T, 



" int 



(33) 



+ 862(p,-b)^(k,-b)^Ti^) 



int ■ 



All e^-terms are at 6 ~ i? of the same order of value, 0(1), and are small compared to the 
main term in Eq.(30), which is of order 0{E?) ~ A^/^. So extra unfactorized correction to 
Eq.(31) due to diagram 7b is inessential. 
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The situation is different for diagram 4c. Its isotropic part has a smallness of order A~^/^ 
in comparison with isotropic part of -f"/^]/- anisotropic part contains a 'non-flow' term 
{pt ■ k^)^ which contains the overlap function Tab (not its derivatives). As a result its contri- 
bution to (31) is of the same order as due to factorized structure (ptb)^(ktb)^ from Eq.(30). 
This contribution violates factorization in Eq.(31) for two-particle correlation. Note that 
the 'flow' and 'non-flow' structures have different dependence on impact parameter. 

6 Role of multipomeron interactions 

In this section we try to estimate contributions to azimuthal asymmetry coming from mul- 
tipomeron interactions. The simplest reggeon diagram of Fig. lb gives main contribution 
to inclusive cross section. Contributions of multipomeron diagrams of Glauber type cancel 
each other because of the AGK rules [9] ®. Note that in 'event-by-event' analysis the number 
of secondary particles is fixed for every individual event, and the AGK rules may not be 
fulfilled. However, all extra interactions change both isotropic and anisotropic parts of the 
amplitude to the same extent, and therefore do not change the flow coefficients Vn- 

In reggeon theory with ap{0) > 1 the number of pomeron exchanges increases with 
energy, so a role of inter-pomeron interactions grows. Detailed analysis of multipomeron 
effects will be published elsewhere, and here we will discuss briefly physical reasons for 
possible increase of azimuthal anisotropy. 

When energy increases a longitudinal length of the nucleon tube with fixed transverse 
position is contracted and multiperipheral fluctuations connected to each of fast nucleons 
overlap each other in longitudinal size. Therefore a probability for particles (partons) of 
different chains (related to different nucleons) to interact increases. Two chains can either 
fuse to a single chain (triple pomeron interaction) or rescatter without changing a whole 
chain number (2-2 or, more generally, n-n pomeron vertices). Effects of fusion become 
noticeable already at RHIC energies, resulting in reduction of secondary density about 
two times as compared with the eikonal approximation [10]. This effect is mainly due 
to diagrams with three-pomeron interaction (Fig. 8a). But due to their sign- alternating 
character a profile of the amplitude becomes more smooth and the azimuthal asymmetry 
becomes even smaller because of smaller gradients of parton density. 

Different situation is realized for diagrams with chain rescattering (Fig. 8b) which give 
positive contributions to inclusive cross section. Even if the total inclusive cross section 
changes slightly due to such rescatterings, its anisotropic part may be strongly increased. 
Indeed, the expression (18) for inclusive cross section will contain in case of n interacting 
chains a convolution T4 C?) instead of Ta®Tb- Azimuthal asymmetry is determined by 
density derivatives (see (21)), therefore an extra factor {n{b))'^ will appear in expression 
for V2- The number of interacting pomerons at RHIC energies is not large ((n(6)) ~ 2 for 
central coUision) , but increase of V2 can be considerable because the effect is quadratic in n 
(see Fig. 9 in section 7). Note that {n{b)) increases with atomic number A and changes the 
A dependence for V2 also. The same mechanism may give an increase of (p^). 

^ The cancellation occurs between diagrams with different numbers of reggeon chains (cut reggeons) 

considered together with sign-alternating absorptive corrections (imcut reggeons). Note that consistent 
account for absorptive corrections is an important feature of the reggeon approach which differs it from all 
classical considerations of multiparticle production. 
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(a) (b) 

Figure 8: Multipomeron diagrams for inclusive production: (a) - diagrams with fusion and 
splitting of parton chains; (b) - diagrams with rescattering of n parton chains. 



7 Numerical estimates 

Let us present examples of numerical estimates of elliptic flow for the exponential parametriza- 
tion (6) of the inclusive vertex. They are illustrative only. We intentionally do not fit exper- 
imental data because, first, little quantitative information on the inclusive vertex 7 exists 
and, second, it is necessary to have more reliable model for nuclear multipomeron effects. 
In this publication we aim to demonstrate that qualitative features of elliptic and directed 
flows peculiar to the mechanism under discussion are in agreement with experimental ob- 
servations. In particular, as it is seen from Eq. (21) the value of V2ipt, b) grows with increase 
of both transverse momentum pt and impact parameter b. 

v2 vs pt at different centrallties 

^ 0.2 
5 

° 0.175 

o 

9. 0.15 
0.125 
0.1 
0.075 
0.05 
0.025 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

tronsverse momentum p, (GeV/c) 

Figure 9: Elliptic flow f 2 as a function of transverse momentum pt in Au — Au collisions for 
two centrality values: C = 13% (6 ~ 5.6 fm) - curve 1; C = 27.5% (6 ~ 7.9 fm) - curve 2. 
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elliptic flow vs centrolity 




Figure 10: Elliptic flow V2 as a function of impact parameter b at pt — 0.5 GeV/c. Calcula- 
tions for 6 > 8 fm are shown with dotted curve. Thick solid line corresponds to calculations 
at e = 0.4. Thin solid line corresponds to calculations at e = 0.1, and dotted line - at e = 0.1 
multiplied to the mean rescattering number squared n(6e//) = n(0.65 6). The dashdotted 
line at this picture gives n(6e//) as a function of b (the right scale). 



It is shown on Fig. 9 pf-dependence of elliptic flow for two values of centrality. Calcu- 
lations were made for two different values of centrality at e = 0.4 for the standard Woods- 
Saxon parameterization of nuclear density. Note that for large pt the growth of V2 is stopped 
due to higher order of e though it is not quite justified to use Eq.(6) in high pt region. Ex- 
perimental data are taken from [11]. The too fast increase of theoretical predictions with 
centrality may be related to a neglect of multipomeron diagram contribution (see below). 

Fig. 10 presents the 6-dependence of eUiptic fiow. According to Eq. (21) V2 grows as 6^ 
(thick solid line, e = 0.4) in contradiction with the experimental behavior at very large b 
when nuclei only touch each other. Let us note here that in the case of very large impact 
parameter values where nuclear densities are quite small an accuracy of the model is low. 
This is why calculation results for 6 > 8 fm are shown with dotted curves. The dashed curve 
corresponds to the simplest account of nuclear effects due to inter-pomeron interactions - 
calculations at smaller parameter value e = 0.1 (thin solid line) are multiplied to the factor 
n(6e//)^- The function n(6e//) gives an effective number of interacting pomeron ffuctuations 
^. It is also presented at the figure (dash-dotted line) and one can see that the maximal 
number of pomeron interaction is about 2 for central collisions. Experimental data are from 
[11]. 

Dependence of the directed fiow vi on rapidity y calculated at 5 = 0.1 for two energies 
(200 and 64 GeV) is presented on Fig. 11. It is antisymmetric and concentrated at the 
edges of the rapidity interval rapidly decreasing to the middle of it. The experimental data 
corresponding to centrality 30 - 60 % are taken from ref.[12]. 

^ The function n{b) was estimated in the eikonal approximation with effective parton-nucleon cross 
section value a ~ 10 mb. The mean position of interacting nucleon in a nucleus was taken as beff = 0.65 b. 
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directed flow vs rapidity 




8 Conclusions 



The main aim of this paper is to draw attention to the mechanism of appearance of azimuthal 
correlations between produced particles related to specific parton correlations in initial states 
of colliding nuclei, and which is not caused by particle interactions in the final state. 

The mechanism under discussion is connected with correlation between parton trans- 
verse momentum and its transverse position inside nucleus. Transverse momenta of partons 
in incoming nuclei (before collision) are already partially aligned along transverse nuclear 
radiuses. The non-central collision of nuclei selects the asymmetric overlapping part of the 
fast nucleus resulting in — b correlations for produced particles. 

We discussed this phenomenon using the simplest reggeon diagram in order to empha- 
size the fact of its existence. For quantitative description of asymmetries it is necessary to 
account for a more complicated mechanism related to interactions between different chains. 
We emphasized that this mechanism can substantially increase elliptic flow and modify 
its atomic number dependence. Role of multipomeron interactions describing these effects 
increases with energy and is important at RHIC energies. This mechanism influences de- 
pendence of asymmetries on centrality, transverse momenta and rapidity. 

Evidently the final state interaction can also contribute to azimuthal asymmetry of 
produced particles. In order to obtain an information on properties of hadronic matter at 
high temperatures and densities it is crucial to separate effects due to initial and final state 
interactions. In particular it is important to carry out a detailed investigation of azimuthal 
correlations of particles produced in nucleon-nucleon and nucleon-nucleus interactions. Such 
studies will be possible in high-multiplicity events at LHC 
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